Astronomy & Astrophysics manuscript no. 1713 


© ESO 2009 


June 22, 2009 





On 
O 
O 
(N 

C 

(N 
(N 



6 

o3 



in 

m 

O 

On 
O 



Theoretical amplitudes and lifetimes of non-radial solar-like 

oscillations in red giants 

M.-A. Dupret 1 ' 2 , K. Belkacem 2 ' 1 , R. Samadi 1 , J. Montalban 2 , 0. Moreira 2 , A. Miglio 2 , M. Godart 2 , P. Ventura 3 , 
H.-G. Ludwig 4 ' 5 , A. Grigahcene 6 , M.-J Goupil 1 , A. Noels 2 , and E. Caffau 5 

1 LESIA, Observatoire de Paris, CNRS UMR 8109, Universite Paris Diderot, 5 place J. Janssen, 92195 Meudon, France 

2 Institut d'Astrophysique et de Geophysique de l'Universite de Liege, Allee du 6 Aoflt 17 - B 4000 Liege, Belgium 
INAF - Osservatorio Astronomico di Roma, MontePorzio Catone (RM), Italy 

4 CIFIST Marie Curie Excellence Team 

5 GEPI, Observatoire de Paris, CNRS, Universite Paris Diderot; 92195 Meudon Cedex, France 

6 Centro de Astrofisica da Universidade do Porto, Rua das Estrelas, 4150-762 Porto, Portugal 



Preprint online version: June 22, 2009 

ABSTRACT 

Context. Solar-like oscillations have been observed in numerous red giants from ground and from space. An important question arises: 
could we expect to detect non-radial modes probing the internal structure of these stars? 

Aims. We investigate under what physical circumstances non-radial modes could be observable in red giants; what would be their 
amplitudes, lifetimes and heights in the power spectrum (PS)? 

Methods. Using a non-radial non-adiabatic pulsation code including a non-local time-dependent treatment of convection, we compute 
the theoretical lifetimes of radial and non-radial modes in several red giant models. Next, using a stochastic excitation model, we 
compute the amplitudes of these modes and their heights in the PS. 

Results. Distinct cases appear. Case A corresponds to subgiants and stars at the bottom of the ascending giant branch. Our results 
show that the lifetimes of the modes are mainly proportional to the inertia /, which is modulated by the mode trapping. The predicted 
amplitudes are lower for non-radial modes. But the height of the peaks in the PS are of the same order for radial and non-radial modes 
as long as they can be resolved. The resulting frequency spectrum is complex. Case B corresponds to intermediate models in the red 
giant branch. In these models, the radiative damping becomes high enough to destroy the non-radial modes trapped in the core. Hence, 
only modes trapped in the envelope have significant heights in the PS and could be observed. The resulting frequency spectrum of 
detectable modes is regular for £ = and 2, but a little more complex for I = 1 modes because of less efficient trapping. Case C 
corresponds to models of even higher luminosity. In these models the radiative damping of non-radial modes is even larger than in the 
previous case and only radial and non-radial modes completely trapped in the envelope could be observed. The frequency pattern is 
very regular for these stars. The comparison between the predictions for radial and non-radial modes is very different if we consider 
the heights in the PS instead of the amplitudes. This is important as the heights (not the amplitudes) are used as detection criterion. 

Key words. Stars: oscillations - Stars: interiors - Convection 



while others are more complex. Modes with very long lifetimes 
(more than 50 days) are detected. 



•i-H 1. Introduction modes with very short lifetimes in f Hydrae. We finally mention 

/\ ' the recent detection by the CoRoT satellite of solar-like radial 

H ' _ .,, . . , . , ... , , , . „ , and non-radial oscillations in many red giants (De Ridde ret al.l 

03 Oscillations of high radial order p-modes have been firmly . , . , „ . c a ■ ^ a » 

, a j • 1 j r j- i i • j 2009). Different kinds of power spectra are found in these data 

detected in several red giants, from radial velocity data ^ — , , .. r „ , . . . 

it; — j ^ i i -./-./-v^ Li r. a a- — \ r . ,. — i -,/-w-yA j . or unprecedented quality, some of them showing regular patterns 

Frandsen et al. 2002; Bedding & Kieldsen 2006) and photomet- ... r n J , , V1 _ r ,- £ ,. 

ric data from space (IBarban et al.ll2.007c De Ridder et al.l l2009). 
These oscillations are stochastically excited by turbulent mo- 
tion in the outer layers of th ese stars. A part i cularly interest- From a theoretical point of vie w, an important study was 
ing case is e Oph, in which iDe Ridder et al.l d2006l) detected made by iDziem bowski et al. (2001) (hereafter D01). They un- 
equidistant frequencies from ground-based spectroscopic ob- derline the mode napping phenomenon present in these stars and 
servations. iHekker et al.l d2006l) analysed the line-profile varia- the role of radiative damping for several models with different 
tions and advocate the presen ce of non-radial modes. For the luminosities. They predict that non-radial modes not trapped in 
same star, Barban et al. (2007) used data obtained by the MOST the en velope would have much smaller amplitudes than radial 
satellite. They detected at least 7 e quidistant peaks and de rived modes. iHoudek & Goughl d2002l) carried out theoretical compu- 
mode lifetimes of about 2.7 days. Kallin ger et al.l d2008al) re- tations of the amplitudes and lifetimes of radial modes in the 
examined the available data for e Oph and claim to detect at star £ Hydrae. They found lifetimes of the order of 15-20 days 
least 21 indep endent modes with lon ger lifetimes between 10 with a bump goi ng up to infinite value (u nstable mode) around 
and 20 days. Kallinger et alj d2008bl) claim to have detected 1 10 (/Hz. Finally. Christ ensen-Dals gaard (2004) also discuss the 
non-radial modes with long lifetimes (more than 10 days) in mode trapping and inertia of non-radial adiabatic modes in red 
the star HD 2 0884 observed by M OST. A very different result giants and subgiants and the possible consequences for the am- 
was found by Stel lo et all Q2006) who claim to have detected plitudes. 
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In the present study, we compute theoretical lifetimes of 
radial and non-radial modes of several red giant models, us- 
ing a fully non-radial non-adiabatic pulsation code includ- 
ing the non-loca l time-d ependent treatment of convection by 
Grig ahcene et al.l (|2005) (G05). We discuss how these life- 
times are affected by the mode inertia, the radiative damping 
in the core and the convective damping in the upper part of 
the convective envelope . Using a stochastic excitation model 
dSamadi & Goupill f2001; Samadietal. 2003b a), we also com- 
pute the amplitudes of radial and non-radial modes and their 
theoretical heights in the power spectrum (PS). As we show be- 
low, predictions for the heights are very different from the ampli- 
tudes. Although non-radial modes in general have smaller am- 
plitudes than radial modes, this is not always the case for the 
heights. As the heights in the PS are used as a criterion for the 
detection of frequencies, it is important to make this distinction. 

In Sect. |2] we discuss the different energetic processes de- 
termining the driving and damping of the modes: the radiative 
damping in the core, the coherent interaction of convection with 
oscillations and the stochastic excitation by turbulent motion in 
the upper part of the convective envelope. In Sect. [3] we present 
the numerical tools, models and method used to solve the prob- 
lem. In Sect.|U we present our results for specific red giant mod- 
els representative of the different cases that can occur. 



2. Energetic aspects of oscillations 

Different mechanisms lead to the driving and damping of oscil- 
lations. For observed modes, two possibilities exist. Firstly, they 
could be self-excited. In such case the amplitudes are expected to 
grow until a large amplitude limit cycle is reached. The observed 
amplitudes of the high radial order p-modes in red giants are of 
several m/s. This is too small for such an interpretation and we 
reject this possibility here. Secondly, the modes can be damped, 
but turbulent motions supply them with energy in a stochastic 
way, allowing them to reach observed amplitudes. We assume 
here that the high frequency modes of red giants are stochastic 
excited modes. The damping rate of the modes is given by an 
integral expression of this type: 



v = - 



2<tI\UR)\ 2 M 



(1) 



where we assumed the time-dependence exp(/<xf - rjt) (t is the 
time, cr is the angular frequency and 77 is the damping rate). J dW 
is the work performed by the gas during one oscillation cycle. £ 
is the displacement vector, / is the dimensionless mode inertia: 







|£| 2 dm/(M|&(*)l 2 ), 



(2) 



M is the total mass and £ r (R) corresponds to the radial displace- 
ment at the layer where the oscillations are measured (Rosseland 
optical depth tr = 0.1 in our results). In radiative zones, the 
work is easily modeled. It just results from the heat given to 
the gas by radiation during each oscillation cycle. In convec- 
tive zones, the problem is much more complex and different 
terms can lead to damping: variations of convective flux, tur- 
bulent pressure, viscosity, and dissipation of turbulent kinetic 
energy (G05). Stochastic excitation by turbulent motion also oc- 
curs, leading to the observed amplitudes. This is discussed in 
Sect.O 



2.1. Radiative clamping 

In the g-mode cavity of red giants, the radial wave-number k r = 
+ Y)NI(crr) becomes huge because of the high density con- 
trast between the core and the envelope (N is the Brunt-Vaisala 
frequency). Large variations of the temperature gradient ensue, 
with loss of heat in the hot phase, leading to radiative damping. 
In the asymptotic limit, a simple expression for this damping 
can be obtained li5zi embowskil ll977bl IVan Hoolst et al. 1998; 
iGodart et alj2009b : 



p dW 



dr 



K [€(€+! )] 

2<T 3 



3/2 



r 



Vad ~ V Vad^VgL 

V p r 5 



dr, (3) 



where ro is the radius of the convective core (if present), r c is 
the upper radius of the g-mode cavity, K is a normalization con- 
stant obtained by appropriate matching with the envelope solu- 
tion, V a d and V are the adiabatic and real gradients, g the gravity, 
L the local luminosity and p the pressure. In the central regions 
of red giants, Ng/r 5 is very high because of the high density con- 
trast. But at the same time, p increases rapidly as we enter in the 
dense pure He core. As a compromise between these two tenden- 
cies, the integrand of Eq. O) reaches its largest values around the 
bottom of the H-burning shell. The main radiative damping thus 
occurs in this region for red giants. 

On the other side, the contribution of the g-mode cavity to 
the denominator of Eq. ([TJ is simply: 



T 

•Jrr, 



\£\ 2 r 2 pdr 



4nK 



f 



k r dr . 



(4) 



which also can be significant because of the large N (k r cc N 
in the core). Eq. © scales like {R/R c ) 1/2 (M/M e ) I/2 , where 
i"kh = GM 2 /(LR) is the Kelvin-Helmholtz time, R c is the ra- 
dius of the He core, R is the total radius, M c is the mass of 
the He core and M is the total mass. Equation scales like 
(M c /M) l/2 (R/R c f /2 . As the star climbs the red giant branch, L 
and R/R c increase so that the numerator increases more quickly 
than the denominator and the radiative damping of non-radial 
modes increases. 



2.2. Coherent interaction with convection 

A very important region for the driving and damping of the 
modes is the transition region where the thermal relaxation time 
is of the same order as the oscillation periods. For solar-like os- 
cillations, this transition region corresponds to the upper part of 
the convective envelope. Moreover, in this region the time-scale 
of most energetic turbulent eddies is of the same order as the os- 
cillation periods. Hence, the changing physical conditions due to 
the oscillations lead to periodic variations of the turbulent quan- 
tities (convective flux, Reynolds stress, . . . ), which contribute to 
the work (positive or negative) performed by the modes. This 
is what we call the coherent interaction between convection and 
oscillations. It plays a major role in the damping of high radial 
order p-mode oscillations of red giants. 

This interaction is difficult to model and only a few theo- 
ries have been proposed for linear oscillations. Two of them 
make use of a Mixing-Length formalism. The first by Gabriel 
( 19960 an d furth er developed by G05 follows the original ideas 
of lUnnol (119671) . where a turbulent viscous term opposite to 
the buoyancy is introduced. This treatment implemented in our 
non-radial non-adiabatic pulsation co de is used in t his study. 
The second theory was developed by iGoughl dl977bl) and fol- 
lows the "kinetic of gas" picture of the MET. Both theories 
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can include a non-local treatment (fSpiegelll 1 963t lGough|[r977at 
lBalmforth|[T992l: iDupretet al.ll2006bl) . A third formulation, no 
longer based on a mixing-len gth approach but on a Reynolds 
stress one was also proposed bv lXiong et al.l(l 19971) . These theo- 
ries encountered some succes ses, for example they obtain the red 
edge of the instabi lity strip (lHoudekl2000tlXiong & Deng|200ll 
iDupret et al.ll2.Q04l) . The rep roduction of the Solar mode lifetime 
is also possible but n ot easy dBalmforthll992tlXiong et al.l2000t 
IDupret et al.ll2006al) . 





M/Mq 


Tefi 

(K) 


log(L/L ) 


logg 


R/Rq 


(MHz) 


A 


2 


5264 


1.32 


3.26 


5.50 


315 


B 


2 


4892 


1.80 


2.65 


11.0 


82.1 


C 


2 


4665 


2.10 


2.27 


17.2 


34.6 


D 


3 


5091 


2.00 


2.70 


12.8 


88.6 


E 


3 


5222 


2.00 


2.74 


12.2 


96.3 



Table 1. Global parameters of our models. 



2.3. Stochastic excitation 

For stochastically excited modes, the local squared amplitude of 
velocity variation at the layer where it is measured (we adopt 
the optical depth t r - 0. 1 in our resu lts) is given by (e.g., 
Baud in et al . 2005; lBelkacem et al.ll2006h : 



n 



2t]MI 2t]MI 2 ' 
P is the power stochastically supplied to the modes: 



P = 



8M 



(5) 



(6) 



where C 2 , and C 2 S are the turbulent Reynolds stress and en- 
tropy contributions, r espectively (see ISamadi & Goupill 1200 it 
Belkace m et al.l I20061 for details). Both play a significant role 
in red giants and are included in our study. 77 and I are the damp- 
ing rate and inertia defined in Eqs. (HJ and (0, respectively. To 
isolate the effect of inertia in Eq. (0, we introduce the product 
Tl - PI which is, according to Eq. (0, independent of inertia. 

We also introduce the maximum height of the mode profile 
in the PS, which is an observable and will permit us to draw 
conclusions about the mode detectability (see Sect.|4]for details). 
To this end, one has to distinguish between two cases, namely 
resolved and unresolved modes. The resolved modes present a 
Lor entzian profile in th e PS and their h eights are given by (see 
e.g. lChaplin et al .1l2005tlBelkacem et al.ll2Q06h : 



H = 



V 2 (R) 



n 



rir 2 



2 j] 2 MI 2 2MI 2 



(7) 



where t = I/77 is the mode lifetime. 

When t > r o b s /2 (T^ being the duration of observations), 
the modes are not resolved. In the li mit t — » 00, the heights 
in the PS tend to b ehave like (see e.g. iBer thomieu et al. 2001; 
lLochard et aT1l2005l) : 



H m — 



^obs 



V 2 (R) = 



nr obs T 

4MI 2 



(8) 



In the theoretical predictions of the heights presented in this pa- 
per, we choose the very favorable case of the C0R0T long runs 
with r o b s = 150 days. We use Eq. when r < r o b s /2 = 75 days 
andEq. (0whenr > r o b s /2 = 75 days, so that H is a continuous 
function of t. 

Finally, the observed velocity amplitudes are obtained by in- 
tegrat ion of the projected l ocal velocity over the visible stellar 
disk (Dziemb owskj|[l9~77al) . This introduces a visibility factor 
depending on the inclination angle of the star (through the factor 
P'"(cos /)). We do not include this factor in our theoretical study 
(the inclination angle is unknown). 
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Fig. 1. HR diagram with the location of our five models. 

3. Numerical tools, models and method 

3. 1 . Structure models 

The equilibrium models analyzed in this study were computed 
with th e version updat ed for asteroseismology of the code 
ATON dMazzitellilll979l) . A detailed description of the numeri- 
cal tec hniques and implemented physics is given in l Ventura et al.l 
(2008) and references therein. For the models in this paper 
we adopted the clas sical mixing-length theory of convection 
(Bohm -Vitenselll958h . with a value of the mixing-length pa- 
rameter ccmlt = 1-8, and we assume core overshooting that 
was treated as a diffusive mixing process. For the parameter 
describing the extension of th at extra mixing we took ffov = 
0.015 (see IVentura et alJ 1 19981 for details), and for the chem- 
ical composition X = 0.72 and Z = 0.012. The radi ative 
opacities are those from OPAL (Iglesias & Rogers! 1 199*61) for 
the metal mixture of iGrevesse & Noelsl (1 19931) completed with 
lAlexander & Fergus on ( 1994 ) at low temper atures, and the con- 
ductive ones come from lPotekhin et al1(|l999[). Thermodynami c 
quantities are derived from OPAL dRogers & Nayfonovl l2002). 



I Saumon et al.1 d 19951) for the pressure ionization regime and 
Stolzmann & Bloecker (1996) treatment is used for He/C/O 
mixtures. The nuclear cross-sec tions are from the NACRE 
compilation ( Ang uloet al.l |T999), and the surface boundary 
conditions ar e prov ided by a grey atmosphere following the 
iHenvev et al.l d 1965b treatment. Turbulent pressure is not in- 
cluded in our structure models. 

We focus our study on five well chosen models representa- 
tive of the different situations that can occur. Some of their global 
characteristics are given in Table Q] v cut = (g/2) ^JT\pJp is the 
cut-off frequency of an isothermal plane-parallel atmosphere (es- 
timated here at the photosphere). Evolutionary tracks with the 
location of these models in the HR diagram are given in Fig.Q] 
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3.2. Non-radial non-adiabatic oscillations 

We use the non-adiabatic pulsation code MAD (Dupret 2002|) to 
compute the theoretical lifetimes and inertia of the modes. We 
use the TDC treatment of G05 for the modelling of the convec- 
tive flux and turbulent pressure variations entering in the non- 
adiabatic equations. Different ingredients must be specified. We 
keep a local treatment for the convective flux. The variations of 
the turbulent pressure and entropy gradient are treated in a non- 
local way, with values of the non-local parameter s b = 3 and 
c = 3.5 according to the definition of Balmforth (1992). Our 
value of b is derived from a fit of the turbulent pressure in the 
subadi abatic atmosp heric layers of a solar hydrodynamic simu- 
lation dDupret et al.ll2006bb . We recall that our structure models 
are built without turbulent pressure. We compute it a posteriori in 
order to be able to determine the turbulent pressure perturbation. 
This implies that the effect of turbulent pressure on the stratifi- 
cation is not taken into account in our models. However, the im- 
pact of this inconsistency on the theoretical lifetimes is expected 
to be negligible compared to others aspects. The main source of 
uncertainties in any TDC treatment comes from the perturbation 
of the closure equations. In our TDC treatment, this uncertainty 
appears in the form of a free parameter /3, as introduced by G05. 
We take by default the value fi = —2 i. Adopting this value for 
the Sun leads to predictions of the mode lifetimes in reasonable 
agreement with observations. It also gives in all our red giant 
models stable high radial order p-modes, in agreement with our 
working hypothesis of stochastic excitation. However, we dis- 
cuss in Sect. l4.6l the impact of changing this value. 

Some changes were required in the MAD code to address 
specific numerical difficulties associated with red giants. First, 
we changed the formalism to a Eulerian one in the core. In 
the g-mode cavity of red giants where <x << Lf, we have 
\p'lp\ « \5plp\ and thus dp/p - (d In pjdf) g r . Hence, dp 
and ^ r are not independent functions and they do not allow us 
to capture the oscillatory behaviour of the eigenfunctions in the 
cavity. Therefore, a Lagrangian formalism lacks numerical pre- 
cision in the dense region of evolved stars, which justifies our 
switch to a Eulerian formalism. 

Another significant numerical difficulty comes from the very 
large number of nodes of th e eigenfunctions in the g-mode cav- 
ity. IVan Hoolst et alj dl998l) and D01 used a mixed treatment 
to address this difficulty: numerical solutions of the fully non- 
adiabatic equations in the envelope are matched with asymptotic 
solutions in the g-mode cavity. Our way to proceed is slightly 
different. Firstly we solve the fully non-adiabatic equations from 
the center to the surface, which gives a first estimation of the 
lifetimes. Second, we separate the integrals in Eqs. (Q]) and (O 
in two parts: region 1 where cr « L(,N and region 2 (the re- 
maining). In region 1 we use the results obtained in the asymp- 
totic quasi-adiabatic limit (Eqs. (f3]l and (HJl), and in region 2 we 
use the fully non-adiabatic eigenfunctions obtained numerically 
in the first step. We then compare the results obtained by the 
two approaches: in many cases they are very close, showing that 
our fully numerical solution of the problem is precise enough. 
In other cases corresponding to high luminosity red giants, the 
number of nodes becomes huge (> 1000), which implies loss 
of numerical precision. In such cases we are more confident in 
using results with the asymptotic treatment in the g-mode cavity 
and the fully non-adiabatic solution in the rest of the star. 

A last important numerical issue is associated with the mode 
trapping phenomenon, which is discussed in detail in Sect. 14.11 
Algorithms solving the non-adiabatic problem converge towards 
the closest eigenvalue in the complex plane (the real part corre- 



sponds to the angular frequency cr and the imaginary part to the 
damping rate rj). Some non-radial modes trapped in the envelope 
have much lower inertia than the others. Hence, their damping 
rates are much higher than for other non-radial modes, higher 
than the frequency separation between consecutive non-radial 
modes: rj > cr n -\ - cr„. Hence, with a real initial guess of the 
eigenvalues (e.g. the adiabatic frequencies), it is impossible to 
converge towards the modes strongly trapped in the envelope. 
The only way is to adopt a complex initial value; for exam- 
ple, interpolating the complex non-adiabatic eigenvalues of ra- 
dial modes gives adequate initial values, allowing the solution to 
converge towards the non-radial modes strongly trapped in the 
envelope. 

3.3. Ingredients for the stochastic excitation models 

The power stochastically supplied to t he modes, P in Eqs. j5] ) 
and (O, is computed as described in Belkacem etalJ (2008). 
The typical convective length-scales as well as the kinetic en- 
ergy spectrum are poorly known for red giant stars. Given 
that the predictions of the stochastic excitation models strongly 
depend on them, we infer both the kinetic energy spectrum 
and the injection length-scale from a representative 3D nu- 
merical simulation. To this end, we computed a 3D radiation- 
hydrodynamical mo d el atmosphere with the code CO 5 BOLD 
dFrevtag et al.l 120021: IWedemever et al.l 120041) . The simulation 
has a gravity logg = 2.5 and an effective temperature of T e g = 
4960 K. The model has a spatial mesh with 160 x 160 x 200 
grid points, and a physical extent of the computational box of 
573 x 573 x 243 Mm 3 . Using this 3D numerical simulation, we 
determine the injecti on length-scale in the layer where the driv- 
ing is the largest (see lSamadi et al.ll2003bl for details). This esti- 
mation of the injection scale gives values larger than the mixing- 
length by a factor of 5-10 depending on the model. Using the 
mixing-length as the injection scale would give theoretical am- 
plitudes much lower than typical observations. We assume here 
that the injection scale is constant in the driving region. In the 
case of the Sun, such approxima tion underestimate s by about 
25 % the theoretical amplitudes (ISamadi et al.l2003bl) . The mod- 
els considered in our study have different gravities. This affects 
the injection length-scale and we assume here that it scales as 
the pressure scale-height at the photospher e. Such a procedur e 
has been verified in the case of a Cen A (Sam adi et al.l 12008). 
This allows us to compute the excitation rates for all the models 
presented in Sect. 13.11 W e use a Lorentzian pr ofile for the eddy 
time-correlation function (Samadi et al.l2003ah . Eventually, con- 
cerning the ^-dependency of the kinetic energy spectrum (k is 
the wavenumber in the Fourrier space of turbulence), we us e the 
Broad Kolmogorov Spectrum (BKS) dMusielak et al.ll 1994 that 
reproduces the spectrum inferred from the 3D simulation. 

4. Results 

4. 1 . Mode trapping in evolved stars 

Mode trapping is a very important phenomenon allowing to un- 
derst and the behaviour of non-radi al modes in evolved stars. D01 
and|Christensen-Dalsgaard (2004) discussed it in the case of red 
giants. Here we recall basic aspects that are important for the in- 
terpretation of our results. The Brunt-Vaisala (N) and Lamb (L() 
frequencies are huge in the dense core of red giants. Because of 
these huge values, all non-radial modes with frequencies lower 
than the cut-off frequency and larger than the fundamental ra- 
dial mode have a mixed character. In the envelope, cr 2 > L 2 ,N 2 



M.-A. Dupret et al.: Theoretical amplitudes and lifetimes of non-radial solar-like oscillations in red giants 



5 




to 
■a 



Fig. 2. Energy density |d//dlogr| for the mode £ — 2, v = 
55.72 //Hz trapped in the envelope (grey) and the mode 1-2, 
v - 53.87 yuHz trapped in the core (black), for model B. 



■g 0.001 




30 35 40 45 50 55 60 65 70 
Frequency (|iHz) 

Fig. 3. Inertia of £ — 0, 1 and 2 modes of model B. 



and they behave like acoustic modes (as in solar-like main se- 
quence stars). But a 2 < L 2 ,N 2 in the core and the same modes 
behave there like gravity modes. The presence of an evanescent 
region between the p- and g-cavities is at the origin of the mode 
trapping. Some modes have significant kinetic energy in the g- 
cavity and low ones in the p-cavity; they are trapped in the core. 
Others have low kinetic energy in the g-cavity but high ones in 
the envelope; they are trapped in the envelope. 

An illustration of the distribution of energy of a mode 
trapped in the core and a mode trapped in the envelope are shown 
in Fig. |2] where we give |d//dlogT| = -|^| 2 dm/(Mdlogr). 
Mode trapping directly affects the mode inertia. We give in Fig. [3] 
the inertia as a function of frequency for £ — - 2 modes. 
We first notice how dense the spectrum of frequencies for non- 
radial modes is, which is a direct consequence of the huge Brunt- 
Vaisala frequency in the core (according to the asymptotic the- 
ory, P nJ - P n -U - 2n 2 /({({ + 1) jN/rdr). Modes trapped in 
the envelope are the local minima with low inertia in Fig. [3] 
The £ — 2 modes trapped in the envelope have lower iner- 
tia than the £ = 1 modes trapped in the envelope, while it is 
the contrary for all other mo des. In agreement with D01 and 
IChristensen-Dalsgaardl d2004l) . the size of the evanescent re- 
gion and thus the efficiency of trapping increases with the de- 
gree £ of the modes. Th is is very different from RR Lyrae stars 
dVan Hoolst et al.| [l998) because there the evanescent region is 




100 150 200 
Frequency (|iHz) 

Fig. 4. Theoretical lifetimes of radial modes of model A. The 
horizontal line corresponds to our choice of the resolution limit 
(75 days). 



larger for £ = 1 in the observed range of frequency. It is very 
interesting to notice that when only the modes trapped in the 
envelope are considered, the frequency spectrum is similar to 
solar-type main sequence stars, with large and small separations 
appearing clearly. The inertia appears directly in Eqs. dT), © 
and (JHJ. Hence, the mode trapping is expected to strongly affect 
the amplitudes, lifetimes, and heights of the modes in the PS of 
red giants. In the next section, we present our main results for 
different models (TableQ] Fig. [T}. 

4.2. Model A: bottom of the red giant branch 

We first consider the lifetimes of radial modes from the fun- 
damental up to the cut-off frequency. We also give the resolu- 
tion limit of 75 days (horizontal line). As can be seen, we can 
separate the modes into two sets. From the fundamental to p7 
(v = 123/^Hz), the lifetimes are very long and such modes could 
not be resolved. On the other side, for larger frequencies the life- 
times of radial modes are below 75 days and these radial modes 
have a chance of being resolved. Current observations of red gi- 
ant solar-like oscillations show th at at least some of the modes 
can be resolved. We note also that Xiong & Deng] d2007h predict 
unstable modes in the lower part of the spectrum. Depending on 
the parameter j3 used in our TDC model, they can be stable or 
unstable. We restrict our study to the second set of modes, in 
which radial modes could be resolved. We propose to call them 
solar-like modes by analogy with the solar case. Such modes are 
detected (as shown in the introduction) and their physical inter- 
pretation as stochastic modes has been firmly established. 

In the top panel of Fig. [5] we give the theoretical lifetimes 
t — 1 / rj obtained for the £ = 0, 1 and 2 modes of model A. 
According to Eq. ([TJ, the mode lifetime r = I/77 oc /. Moreover, 
we see from Eq. © that the heights in the PS are proportional 
to (t//) 2 for resolved modes. This leads us to consider in the 
second panel the mode lifetimes divided by the inertia (t/7). As 
we showed in Fig. [3] the mode trapping leads to an oscillatory 
behaviour of the inertia with the large frequency separation as 
periodicity. The oscillatory behaviour of the lifetimes is a direct 
consequence of the oscillations of I. But when the ratio between 
the two is considered, the oscillatory behaviour disappears and 
the same results are found for radial and non-radial modes. 

In order to interpret this result, we give in Fig. |6]the cumu- 
lated work integral W(m) - ["* dW/(2crIM) for a typical 1-2 
mode trapped in the core (v = 201.3 fiHz). The surface value 
is -rj R 3 ^ 2 /(GM) 1 ' 2 . We see that significant driving and damp- 
ing only occur in the upper part of the convective envelope, but 
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Fig. 5. Theoretical predictions for ( — 0, 1 and 2 modes of 
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responds to our choice of the resolution limit (75 days). Panel 2: 
lifetimes over inertia (t/7). Panel 3: amplitudes (m/s). Panel 4: 
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Fig. 6. Work integral for the mode i — 2, v = 201 .3 fjHz trapped 
in the core of model A. 



there is no significant radiative damping in the core. As the ra- 
dial component of the displacement dominates in the envelope, 
the numerator of Eq. (HJ is the same at a given frequency, what- 
ever the value of t. This explains why r/7 does not depend on 
1. 

In the third panel of Fig.|5J we give the theoretical amplitudes 
of the modes. The amplitudes of radial modes are higher because 
they have less inertia. The oscillatory behaviour of the inertia 
due to mode trapping leads to oscillations of the amplitudes of 
non-radial modes. 

In the bottom panel of Fig. |5J we give the theoretical heights 
in the PS. We put the resolution limit at 75 days (horizontal line 
in top panel). If the mode lifetime value is larger than this value, 
Eq. ([8]) is used instead of Eq. (0 to obtain the height. For re- 
solved modes (H cc (r//) 2 ), the same heights are found for ra- 
dial and non-radial modes, except for the visibility factor (inte- 
gration over the stellar disk), which is not included in our study. 
The unresolved non-radial modes with long lifetimes (mainly 
I — 2 modes trapped in the core) have smaller heights than the 
closest radial mode. We emphasize that the predictions for the 
heights are different from the amplitudes. Although non-radial 
modes have smaller amplitudes than radial modes, the heights in 
the PS are often of the same order. Contrary to the amplitudes, 
the heights are observation dependent. If r o b s is smaller, more 
non-radial modes are unresolved and their heights in the PS are 
smaller. 

Thus, many non-radial modes would be detectable in the PS 
of stars like model A. Moreover, the interaction between the p- 
and g-mode cavities leads to many avoided crossings. Hence, 
the resulting frequency spectrum is very complex, particularly 
for € — 1 modes and its seismic interpretation is not easy. 

4.3. Model B: intermediate in the red giant branch 

We now present the results obtained for model B that is slightly 
higher in the red giant branch (Fig.[T] TableQ]). In the top panel of 
Fig. [7] we give the theoretical lifetimes obtained for the t — 0, 
1 and 2 modes. Again we see the oscillatory behaviour due to 
mode trapping. But considering now the ratio between the life- 
time and the inertia given in the 2nd panel from the top, we find 
very different results compared to model A. Now, oscillatory be- 
haviour of this ratio is present for non-radial modes, particularly 
at low frequencies. To understand this, we compare in Fig. [8] 
the work integral for two 6 = 2 modes, one trapped in the en- 
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Fig. 7. Theoretical predictions for t — - 2 modes of model B. 
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Fig. 8. Work integrals for the mode £ = 2, v — 55.72 //Hz trapped 
in the envelope (top panel) and the mode t = 2, v = 53.87 fjHz 
trapped in the core (bottom panel) of model B. 



velope (v = 55.72 //Hz, top panel) and the other trapped in the 
core (v = 53.87 //Hz, bottom panel). For the mode trapped in the 
envelope, the driving and damping occur in the upper part of the 
convective envelope, as in model A. But for the mode trapped in 
the core, significant radiative damping occurs around the bottom 
of the H-burning shell (log T 7.5). 

In the third panel of Fig. |7l we give the theoretical ampli- 
tudes of radial and non-radial modes of model B. The ampli- 
tudes of non-radial modes trapped in the core are very small be- 
cause of large inertia and radiative damping. The amplitudes of 
non-radial modes trapped in the envelope (local maxima in the 
figure) are smaller than radial modes but not negligible because 
of similar inertia and negligible radiative damping. 

In the bottom panel of Fig. [7] we give the heights of the 
modes in the PS. Non-radial modes trapped in the core have neg- 
ligible heights and would not be detected. In contrast, non-radial 
modes trapped in the envelope have similar heights compared to 
radial modes. 

Thus, only radial modes and non-radial modes trapped in the 
envelope could be observed in the PS of stars like model B. The 
groups of non-radial modes trapped in the envelope are more 
or less separated from each other by a constant large separa- 
tion. Seismic interpretation of the frequency spectrum of such 
stars would thus be easier than in model A. We notice however 
that the trapping is not perfect for t — 1 modes. Hence, a small 
group of I — 1 modes is detectable around each local maximum. 
This would have to be taken into account in any seismic study. 
Considering observed power spectra, this also means that these 
groups of modes would have to be not confused with single large 
line-widths modes. 



4.4. Model C: high in the red giant branch 

We consider now the results obtained for model C, which is even 
higher in the red giant branch than the previous ones (Fig. [TJ. In 
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Fig. 9. Theoretical predictions for I — - 2 modes of model C. 
Top panel: theoretical lifetimes, panel 2: amplitudes (m/s), 
panel 3: heights in the PS. 



the top panel of Fig. [9] we give the theoretical lifetimes. They are 
very different from those of model B. The oscillatory behaviour 
is no longer present for most t — 2 and low frequency i — 1 
modes. As can be seen in the figure, the transitions to envelope- 
trapped modes are very sharp; these modes, known as Strongly 
Trapped in the Envelope (STE), are nearly perfectly reflected at 
the internal turning point (near the top of the H-burning shell). 
We give in Fig.[l0]the work integral for an t — 2 mode with fre- 
quency very close to but different from an STE mode. The very 
high radiative damping illustrated in this figure occurs for all 
non-radial modes of this model, except the STE mode. Because 
of this and the large inertia, the g-mode cavity dominates the in- 



Fig. 10. Work integral for the mode I = 2, v = 22.83 //Hz of 
model C (this mode is near a mode trapped in the envelope). 

tegrals in the numerator and the denominator of Eq. (Q~|) for all 
non-STE modes. Hence Eq. ([3]) and Eq. (0 are good approxi- 
mations for the full integrals and the K constant simplifies when 
the ratio between the two is considered. This explains why the 
lifetimes do not show oscillations in the top panel. In the sec- 
ond panel of Fig. [9] we give the theoretical amplitudes. They 
are very small for all non-radial non-STE modes, because of 
the high radiative damping and inertia. Finally, we give in the 
bottom panel the theoretical heights in the PS. All non-radial 
non-STE modes have negligible heights because of the large ra- 
diative damping. Thus, only radial modes and non-radial modes 
completely trapped in the envelope could be detected in stars like 
model C. 

4.5. Models D and E: H-shell versus He-core burning 

We now compare the theoretical lifetimes, amplitudes and 
heights of the two models D and E with the same luminosity, 
one before (model D) and the other during (model E) the core 
helium burning phase. 

In Fig.[TT] we give the theoretical lifetimes obtained for the 
I — - 2 modes of these two models. The behaviour is similar 
for the two models. 

In Fig. [12] we give the theoretical heights in the PS for these 
two models. No significant difference between the two is found 
except that the heights in the PS are systematically larger for the 
He burning model, because of its higher effective temperature 
leading to stronger stochastic driving. 

4.6. Uncertainties of the non-adiabatic treatment in the 
convective envelope 

As discussed in Sects. l2~2l and 13.21 the time-dependent interac- 
tion between convection and oscillations is subject to large un- 
certainties. In the TDC treatment of G05, a free complex param- 
eter /3 was introduced in the perturbed closure equations. In the 
previous computations, we adopted f3 = -2 i. We now redo the 
computations with /3 = - 1 .75 i. Larger changes of /3 would give 
unstable modes in the considered range of frequencies, which is 
not realistic in view of the typical observed amplitudes of red 
giant high order p-modes. 

With p = -1.75 i, the convective damping in the upper part 
of the convective envelope is larger. As shown in Fig. [13] this 
results in smaller lifetimes and amplitudes of the modes trapped 
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Fig. 11. Theoretical lifetimes for t — - 2 modes of model D 
(pre-He burning, top panel) and model E (He burning, bottom 
panel). 

in the envelope compared to Fig. [7] In the upper part of the en- 
velope, the eigenfunctions do not depend significantly on the de- 
gree i of the modes. Hence the uncertainties associated with the 
TDC treatment affect radial and non-radial modes trapped in the 
envelope in the same way. 

4.7. Non-linear effects? 

Our study is linear, raising the question of whether such an ap- 
proximation is valid here. Considering first the amplitudes of the 
displacement at the photosphere, we find typical values about 
200 times smaller than the pressure scale height for detectable 
modes. Hence non-linear effects in the superficial layers are 
expected to be small. We consider now the oscillations in the 
g-mode cavity. We can confidently neglect non-linear terms if 
\v r dv/dr\ « \d\/dt\. In a cavity with short wavelength oscil- 
lations, this condition is equivalent to \g r \ « l/k r . In the g- 
mode cavity of envelope trapped modes, |£ r | is typically about 
10 4 times smaller than 1 jk r . Hence the risk of wavebreaking due 
to non-linear effects is negligible in the core. 

5. Conclusions 

We have determined the theoretical lifetimes, amplitudes and 
heights in the PS of radial and non-radial modes of several red 
giant models. The predictions appear to be very different de- 
pending on the evolutionary status of the star. What mainly mat- 



Fig. 12. Theoretical heights in the PS for I — 0, 2 modes of 
model D (top panel) and model E (bottom panel). 

ters is the density contrast between the He core and the en- 
velope. During evolution, it increases considerably as the He 
core contracts and the envelope expands. Moreover the Kelvin- 
Helmholtz time decreases because of the increasing luminosity. 
The consequence of these two effects is an increasing importance 
of the radiative damping of non-radial modes as the luminosity 
of the star increases. We can encounter very different cases de- 
pending on the evolutionary status of the star: 

- Case A: Subgiants and/or low luminosity red giants. 

Our model representative of the class is model A at the bot- 
tom of the giant branch for 2Mq. In these stars, the density 
contrast is not large enough to imply significant radiative 
damping of non-radial solar-type p-modes (R c /R - 0.0075 
in model A). The amplitudes of non-radial modes are smaller 
than radial modes because of larger inertia. However, if we 
consider instead the heights in the PS, they are similar for 
radial and non-radial modes as long as they can be resolved. 
As the heights are used as detection criteria in a Fourier 
analysis, many non-radial modes could be detected in the 
PS of such stars. The interaction between the p- and g-mode 
cavities leads to many avoided crossings. The resulting fre- 
quency spectrum is thus very complex. Asteroseismology of 
such stars would be challenging. 

- Case B: Intermediate models in the red giant branch. 

Our models representative of this class are model B, D and 
E with greater luminosity compared to the previous case 
(log(L/Lo) - 1.8 — 2). In these stars the radius ratio be- 
tween the core and the envelope is smaller {R c /R = 0.003 in 
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model B). As a consequence, the radiative damping around 
the bottom of the H-burning shell becomes large enough 
to destroy the non-radial modes trapped in the core. Only 
the non-radial modes trapped in the envelope and the radial 
modes can survive and have similar heights in the PS. The 
non-radial modes trapped in the envelope are more or less 
separated from each other by a constant large separation, al- 
lowing us to build echelle diagrams. However, the trapping 
is not perfect for £ = 1 modes. Hence, a small group of I — 1 
modes is detectable around each local maximum in the PS 
corresponding to an envelope trapped mode. Such structures 
would not have to be interpreted as single large line-width 
modes. Seismic interpretation of the frequency spectrum of 
such stars would be easier than in the previous case. 
Case C: High luminosity red giants. 

Our model representative of this class is model C with 
M - 2Mq, log(L/L0) = 2.1. This case is not very different 
from the previous one, except that the density contrast be- 
tween the core and the envelope is very large (R c /R = 0.002 
in model C). As a consequence, the radiative damping is 
stronger in the core and destroys all non-radial modes except 
those strongly trapped in the envelope. These modes and the 
radial modes are the only ones that could be detected in these 
stars. The trapping is efficient enough in I = 1 modes and 
only single envelope trapped modes are detectable. Seismic 



interpretation of the frequency spectrum of such stars would 
also be easier. 

We also notice that the duration of observation is very im- 
portant, as it determines the frontier between resolved and unre- 
solved modes. Non-radial modes have always much longer life- 
times (r) than radial modes because of their greater inertia. If 
t S T \, s /2, the modes cannot be resolved and have much smaller 
heights in the PS. With too short observations, as is typically the 
case with ground based observations, only radial and strongly 
trapped non-radial modes have a chance of being detected. On 
the other hand, future space missions will provide longer obser- 
vations: e.g. 5 years for Kepler, 3 years for Plato, increasing the 
ability to detect more non-radial modes. 

Asteroseismology of red giants is a very promising field as 
the theoretical predictions depend strongly on the physical char- 
acteristics of both their deep and superficial layers. 
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